Received 2009 December 24 

Preprint typeset using I^T^X style cmulateapj v. 2/16/10 



EVIDENCE OF POSSIBLE SPIN-ORBIT MISALIGNMENT ALONG THE LINE OF SIGHT IN TRANSITING 

EXOPLANET SYSTEMS 

Kevin C. Schlaufman 1 

Astronomy and Astrophysics Department, University of California, Santa Cruz, CA 95064 

Received 2009 December 24 

ABSTRACT 

There has lately been intense interest in the degree of alignment between the orbits of transiting 
exoplanets and the spin of their hosts stars. Indeed, of 26 transiting exoplanet systems with mea- 
surements of the Rossiter-McLaughlin (RM) effect, eight have now been found to be significantly 
spin-orbit misaligned in the plane of the sky. Unfortunately, the RM effect only measures the angle 
between the orbit of a transiting exoplanet and the spin of its host star projected in the plane of sky, 
leaving unconstrained the compliment misalignment angle between the orbit of the planet and the 
spin of its host star along the line of sight. I use a simple model of stellar rotation benchmarked with 
observational data to statistically identify ten exoplanet systems from a sample of 75 for which there 
is likely a significant degree of misalignment along the line of sight between the orbit of the planet and 
the spin of its host star. I find that HAT-P-7, HAT-P-14, HAT-P-16, HD 17156, Kepler-5, Kepler-7, 
TrES-4, WASP-1, WASP-12, and WASP-14 are likely spin-orbit misaligned along the line of sight. All 
ten systems have host stellar masses in the range 1.2 Mq < M* < 1.5 M 0I and the probability of 
this occurrence by chance is less than one in ten thousand. In addition, the planets in the candidate 
misaligned systems are preferentially massive and eccentric. I combine the distribution of misalign- 
ment in the plane of the sky from RM measurements with the distribution of misalignment along 
the line of sight from this analysis, and the result suggests that transiting exoplanets are more likely 
to be aligned along the line of sight than expected given predictions for the spin-orbit misalignment 
distribution of a population produced entirely by planet-planet scattering or Kozai cycles and tidal 
friction. For that reason, there are likely two populations of close-in exoplanet systems: a population 
of aligned systems and a population of apparently misaligned systems in which the processes that 
lead to misalignment or to the survival of misaligned systems operate more efficiently in systems with 
massive stars and planets. 

Subject headings: planetary systems — planetary systems: formation — stars: rotation — stars: 
statistics 



1. INTRODUCTION 

The degree of alignment ip between the orbit of a tran- 
siting exoplanet and the spin of its host star can provide 
insight into the system's formation. In the Solar Sys- 
tem, the angle tp between the planets' o rbits and the 
spin of the Sun is small, about 7° (e.g. Be ck fc Giles! 
I20051) . That fact provided a key piece of evidence 
supporting the earliest models of planet formation in 
the same disk in which the Sun formed. Simple mod- 
els of close-in exoplanet formation involving assembly 
near the ice line and quiescent migration through a 
gaseous protoplanetary disk predict that orbits should 
have only negl igible misalignments with the spin of the 
host star (e.g. iLin et all [19961 llda fc Linll2008a llbT). Al- 
ternatively, simple models without migration but involv- 
ing dynamical scattering between newly formed plan- 
ets and planetesimals at several AU predict that non- 
negligible misalignments ip can occur (e.g. iFord et all 



negligible misalignments (|Thommes fc Lissauerl [2003). 
Moreover, any multiple systems that migrate will likely 
be dynamically stable before the dissipation of their 
gaseous parent protoplan etary disk, but dynamically 
unstable after dissip ation (|Terquem fc Papaloizoull2007l : 
lOgihara fc Ida 2009J). The Kozai effect (IKozailll962D and 



tidal processes might also contribute (e.g. IWinn et all 



2005t iFabrvckv fc Tremainell2007l: iNagasawa et al.ll2008l: 



200lilYu fc Trernaindl200ltlPapaloizou fc Terquemll2001l: 
Terauem fc Papaloizoull2002t iMarzari fc Weidensch illina 
2002tlChatteriee et al.ll2008tljuric fc Tremainell2008f ). In 



reality, there will be so me degree of migration in all mas- 
sive planetary systems (Ward 1997|), and in certain situa- 
tions migration in multiple systems can still lead to non- 
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Pont I2009D . Despite these complications, the frequency 
of misaligned systems will provide another constraint on 
models of planet formation. 

I describe the geometry of tr ansiting exoplanet systems 
according to the conventions of lFabrvckv fc Winnl (j2009f ) 
in Figured] An exoplanet 's orbit is defined by the vector 
rip and is inclined relative to the line of sight at an an- 
gle i p that can be measured through transit photometry. 
Note that the fact that an exoplanet is observed to transit 
requires that i p w 90°. The host star's spin is defined by 
the vector n s and is inclined at an angle i s relative to the 
line of sight. The angle A is the angular distance between 
the projections of n p and n s in the plane of the sky, and is 
observable thr ough the Ross iter- McLaughlin (RM) effect 
(|Rossiterlll924l : IMcLaughlinl[T92l [?)ntaet^^00a . The 
RM effect - the detection of an anomalous blueshift and 
rcdshift superimposed on the stellar reflex in radial veloc- 
ity measurements during transit - has lately been used to 
measure A for several transiting exoplanet systems (e.g. 
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iFabrvckv fc Winn1l2009l and references therein) . Indeed, 
the latest RM measurements described in Triaud et al. 
(2010, submitted) show that eight of the 26 systems with 
RM measurements are significantly spin-orbit misaligned 
in the plane of the sky. Though i s is not constrained by 
the RM effect, if all of i p , i s , and A are known, the calcu- 
lation of the true deprojected angle tp = arccos(n p • n s ) 
between n p and n s is trivial. 

The inclination i s to the line of sight of stellar spin 
is difficult to measure for individual stars. If the ro- 
tation period of a star is known through photometric 
observations of starspots moving across its surface, and 
if the projected rotational velocity v sini has been mea - 
sured, then i s is easy to compute (e.g. lWinn et al.ll2007l ). 
Astrosiesmology meas urements of i« are also possible in 
special circumstances ([Gizon fe Solan ki 2003). Unfortu- 
nately, period measurements are resource-intensive and 
require well-sampled and precise photometric time se- 
ries. They may also be biased toward stars with many 
starspots and therefore biased toward more active and 
younger stars. On the other hand, the projection of stel- 
lar rotation velocity along the line of sight v sin i can be 
measured cheaply by examining the rotational broaden- 
ing of absorption lines in a single high-resolution stellar 
spectrum. If the radius of the star is known from other 
means, it is easy to forward-model the observable quan- 
tity v sin i from the stellar rotation period P* predicted 
from theoretical or empirical models given a distribu- 
tion for i, Monte Carlo simulations can then derive the 
distribution of possible vsmi measurements. That dis- 
tribution can then be compared to observational data for 
transiting exoplanet systems to determine if the obser- 
vations are consistent with the alignment of a transiting 
planet's orbit with the spin of its host star. 

Calculating the rotation period of Sun-like stars as 
a function of mass and time is a tractable problem 
for stars older than about 650 Myr. First, the ro- 
tation period of isolated Sun-like sta rs monotonically 
slows down as they age (jKraftJ I1967D . Second, even 
though a population of Sun-like st ars is formed with 
a broad distribution of p eriods P* (jAttridge fe HerbsD 
H992tlChoi fc Herbstlll996l ). by the time that population 
is the age of the Hyades and Praesepe (about 650 Myr) 
the broad distrib ution has conyerged to a well-define d 
function of mass (iTassoull 12005 llrwin fe Bou^ien 2009). 
Two classes of physical models have been proposed to 
explain this early-time behavior : magnetic break i ng by 
protostellar disks (e.g. IKoenigll 11991b IShu et all Il994t 



iCollier Cameron and accretion-driven stellar 

winds (e.g. iMatt fc Pudritzl I2005D . Though the physi- 
cal process that produces this effec t in the first ~ 650 
Myr is not yet known, the data (IRadick et al.l Il987t 
iProsser et all 119951: iScholz fc Eisloffell 12007ft and t heo- 
retical models (jKeppens et al.1 119951 : iTassoull l2000f ) are 
in agreement on the rotational evolution of Sun-like 
stars after this time. Indeed, Sun-like stars likely lose 
angular momentum through magnetized stellar winds 
(|Weber fc Davislll967t iMestel 119681: lKawaler|[T98l . and 
their rotation slows as P» cx f 1 / 2 , as described observa- 
tionally as wsini oc t~ x / 2 by iSkumanichl (|1972f) . 

In this paper, I use a simple empirical model to cal- 
culate the rotation period of Sun-like stars as a func- 
tion of mass and age. I combine that model with a 



Monte Carlo simulation to marginalize over the uncer- 
tain masses, radii, and ages of a control sample of 866 
stars from the Spec troscopic Properti e s of C ool Stars 
(SPOCS) catalog of lValenti fc Fischer! ([2f)05h to deter- 
mine the expected range of projected rotation vsmi for 
each star in the control sample under the assumption 
that the inclination distribution of the SPOCS sample is 
isotropic. Likewise, I perform a similar calculation for 
the host stars of transiting exoplanets, this time assum- 
ing alignment between the host star's spin and the orbit 
of its exoplanet. I identify those transiting exoplanet 
systems in which the disagreement between the v sin i 
predicted from the simple empirical model under the as- 
sumption of spin-orbit alignment and the observed v sin i 
is larger than the equivalent disagreement for any of the 
866 stars in the SPOCS control sample. In particular, 
those systems that have anomalously small v sin i mea- 
surements relative to the measurement errors and the 
width of simulated v sin i distribution are possibly spin- 
orbit misaligned along the line of sight. I describe the 
details of my calculation in §2. I discuss my results in 
the context of close-in planet formation in §3. I summa- 
rize my findings in §4. 

2. ANALYSIS 

I combine the well defined mass-rot ation period rela- 
tion in the Hyades and Praesepe (e.g. llrwin fc Bouvierl 
I 2009T) with the relation P oc t 1 / 2 (e.g. iWeber fc Davisl 
119671 iMestell 119681 : iKawalerl 119881 ) to compute the ex- 
pected rotation periods of Sun-like stars as a function 
of mass and age. I use a Monte Carlo simulation to 
transform that function into a vsmi S i m distribution by 
marginalizing over the uncertain distributions of stellar 
mass M*, radius P*, age r*, and the applicable incli- 
nation distribution for the sample. As a result, I can 
compare the observed i>sinj & s measurements with the 
mean v sin i 3 i m estimate from the Monte Carlo simula- 
tion to determine the degree of agreement between the 
two. In order to determine the degree of disagreement 
that can be attributed to random effects, I compare the 
measured values of v sin i f, s of Sun-like stars in the Solar 
Neighborhood from the SPOCS catalog with the mean 
v sin i s i m from the Monte Carlo simulation given the stel- 
lar parameters of those stars and the standard isotropic 
inclination distribution. In this way, I use the control 
sample of SPOCS stars to determine a threshold level of 
disagreement that can be expected between the simple 
empirical model and observation given the observational 
uncertainties and the imperfections of the model. I then 
do a similar calculation for the host stars of transiting ex- 
oplanets, this time assuming alignment between the spin 
of the host star and the orbit of its exoplanet. I iden- 
tify those exoplanet host stars in which the predicted 
vsmi S i m disagrees with the observed value i>sinz fc s rel- 
ative to the measurement errors by an amount greater 
then an equivalent disagreement for any of the 866 stars 
in the control sample. Systems with v sin i i, s slower than 
v sin igim by an amount large relative to the quoted er- 
rors in t>sini f, s and the width of the simulated vsini sim 
distribution given the uncertainty in stellar mass, radius, 
and age are either systems with rotation properties unlike 
any star in the control sample or in which the inclination 
of the transiting planet's orbit i p and the inclination of 
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the host star's spin i a are misaligned. Since the SPOCS 
sample is representative of the sample of transiting ex- 
oplanet host stars, the former is very unlikely. For that 
reason, the best explanation of the anomalously slow pro- 
jected rotation is spin-orbit misalignment along the line 
of sight. 

2.1. SPOCS Control Sample 

I exami ne a control samp l e of st ars from the SPOCS 
catalog of Valcnti fe Fischen (|2005f) to quantify the scat- 
ter about the simple mode l. The SPOCS sam ple selec- 
tion is described in detail in lMarcv et al.l (|2004D . In sum- 
mary, the bulk of the sample was initially selected from 
the Hipparcos catalog such that each star has B — V > 
0.55, is no more than three magnitudes above the main 
sequence, is not a known spectroscopic binary, and has 
no stellar companion within two arcseconds. Subsequent 
spectroscopic measurements of Call H and K were then 
used to downselect to those stars with ages r* > 2 Gyr, 
though an additional 100 stars with ages between 50 and 
500 Myr were later added in. As a result, the SPOCS 
sample has both young and evolved s tars as indicated in 
Figure 14 of iValenti fc Fischerl (|2005l ). 

In addition, previously unknown spectroscopic binaries 
were removed from the SPOCS candidate list, so there 
are no spectroscopic binaries in the final catalog. Short- 
period low-mass binary companions or exoplanets would 
have been readily apparent in the original California- 
Carnegie planet search and would have been announced 
as such long ago. For those reasons, the rotation prop- 
erties of stars in SPOCS sample should be unaffected by 
close companions. Moreover, any overlap between the 
SPOCS sample and the hosts stars of transiting exoplan- 
ets is negligible. Collectively, all of these selections and 
observations suggest that the SPOCS sample is a fair 
control sample for this analysis. 

2.2. Detailed Description of the Monte Carlo 

I model the initial mass-period relation of a 650 
Myr population of Sun-like stars b y binning the avail- 
able Hyades and Praesepe data (iRadick et al.l Il987t 
iProsser et aLI [19951: iScholz fc E isloffcl 2007) as summa- 
rized in llrwin fc Bouvierl ([2009) in 0.1 M© bins. The 
llrwin fc Bouvieif ([2009) data only extends up to about 
1.2 Mq, so I supplement it with average rota tion peri- 
ods fo r more massive field stars as presented in IMcNalivl 
(fl965T) . I use natural cubic splines to interpolate the 
binned and supplemented data and find that it minimizes 
the sum of square residuals relative to high-order poly- 
nomial interpolation. I determine the expected rotation 
period of a Sun-like star as a function of mass and age by 
evolving the initial conditi on set by the stellar mass ac- 
cording to the relati on Ce.g. lWeber fc Davislll967tlMestell 
[19681 IKawalen[l988l) 

n(M„r,)=P„ (M,)(^^) 1/2 (1) 

where P*,o(M*) can be approximated as a fifth- 
order polynomial in with coefficients in in- 
creasing order (5.558509,-77.843724,357.387761,- 
502.288156,282.908686,-56.329952). I summarize 
my technique in Figure [2] 



The model described above is likely too simplistic, the 
relevant input parameters (stellar mass, radius, and age) 
are all uncertain, and the output (stellar period) is diffi- 
cult to observe. For those reasons, I use a Monte Carlo 
simulation to determine the range of observable v sin i S i m 
expected for a sample of stars given mass, radius, and 
age as well as the uncertainties in each of those quanti- 
ties plus the appropriate inclination distribution for that 
sample. In order to determine the degree of disagree- 
ment between the v sin i S i m predicted from the simple 
empirical model and the observed v sin z t> s that can be 
expected given the observational uncertainties and the 
imperfections of the simple empirical model, I use the 
control sample of stars from the SPOCS catalog. 

For each star in the SPOCS catalog, I use the esti- 
mated stellar mass and uncertainty M* and o~m (column 
9 of Table 9) , estimated stellar radius and uncertainty i?* 
and an (column 8 of Table 9), and age range At (col- 
umn 14 of Table 9). In the Monte Carlo, I sample each 
star's mass from a uniform distribution in mass between 
M* — om /2 and + <tm/2, its radius from a uniform 
distribution in radius between i?„ — ctr/2 and J?* + <7r/2, 
and its age from a uniform distribution in age between 
the given lower and upper age limits. I sample the in- 
clination to the line of sight of each star's rotation axis 
from the standard random distribution i ~ arccos(l — U), 
where U is drawn uniformly from the interval [0,1] . Using 
these parameters, I evolve the initial period given by the 
randomly selected mass and the initial period-mass rela- 
tion described in Figure [2] to the randomly selected age 
of the system according to Equation [TJ I then compute 
v sin i S i m using the randomly generated radius of the star 
and its randomly selected inclination. I repeat this pro- 
cess 1000 times and thereby derive for each star the dis- 
tribution of possible v sin i S i m values, its mean vsmi S i m , 
as well as the width of the distribution a S i m . I compare 
the predicted v sin i S i m with the observed v sin i Q b s rela- 
tive to the width of the Monte Carlo distribution o~ S i m 
and the error in the observed v sin i b s measurement a b s 
through the rotation statistic 0: 



_ v sin i slm - v sin i obs 

r~2 — ; — 2 — \ > 

V a obs + a sim 

Note that large positive (negative) values of O indicate 
slower (faster) than expected wsini values. Stars that 
are well fit by the model will have small absolute values 
of O. 

I repeat this calculation on each of the 866 stars from 
the SPOCS catalog with M* < 1.5 Mq, assuming the 
standard isotropic distribution of inclination. This dis- 
tribution of O is the null hypothesis that I will use to 
compare to the distribution of O measured in the transit- 
ing exoplanet systems under the assumption of spin-orbit 
alignment. 

I collected the properties of the 80 transiting exoplan- 
ets systems known as of May 2010, of which the 75 listed 
in Table [T] have measured mass, radius, i>sini £, s , and 
orbital inclination. In this case, the stellar masses and 
radii are more precisely measured than those available 
in the SPOCS catalog, and the inclination of each exo- 
planet 's orbit is known from transit photometry. I again 
use a Monte Carlo simulation to determine the range of 
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observable v sin i S i m . I use the measured values for each 
star's mass and radius, and I sample each star's age uni- 
formly between lower and upper bounds for its age from 
the literature. If there are no age estimates available in 
the literature, I sample the age of the star uniformly be- 
tween 650 Myr (the age of the Hyades and Praesepe) and 
its main sequence lifetime. If there are only upper lim- 
its on vsiriiobs in the literature, I assume that the true 
value corresponds to the published upper-limit, as this is 
the most conservative way to approach a search for stars 
with anomalously slow projected rotation. In those situ- 
ations, I assume that the error in vsm.i i, s is 1 km s _1 . I 
assume that the inclination of the host star's spin relative 
to the line of sight i s is the same the inclination of the 
exoplanet's orbit i p . I repeat this process 1000 times and 
compute the rotation statistic O for each stellar host; I 
report the results of this calculation in Table [1] 

2.3. Results 

I plot the distribution of 9 for both samples in Fig- 
ure [3] The O distribution of the control sample peaks 
near = 0, indicating that the simple empirical model is 
accurate in that the expected residual between the model 
and the observations is close to zero. The distribution of 
falls off quickly as increases, and there are no control 
stars from the SPOCS catalog with > 2.9. However, 
there does exist a tail of stars with much larger than 
expected vsini values, and therefore large negative 
values. In other words, some stars in the control sample 
of SPOCS stars are fast rotators; this is to be expected 
in any sample of Sun-like field stars. On the other hand, 
there is no tail of SPOCS stars rotating arbitrarily more 
slowly than the simple empirical model predicts. For that 
reason, host stars of exoplanets that are observed to have 
> 2.9 are apparently unlike any of the 866 stars in the 
control sample - either they rotate more slowly relative 
to the model than any star in the control sample or their 
anomalously slow projected rotation is due to a pole-on 
view of the star and therefore spin-orbit misalignment. 

I define fast rotators as those stars with < —5.2, 
corresponding to the fastest 5% of the distribution for 
the SPOCS control sample. Note that the fraction of fast 
rotators in each sample is statistically indistinguishable - 
about 5% ±3% in the SPOCS sample and about 7% ±3% 
in the transiting exoplanet host star sample. The consis- 
tency of the two samples in the fraction of stars that are 
fast rotators is in sharp contrast to the disparity at the 
other end of the distribution. Also, the peak of the 
distribution of transiting exoplanet systems is offset to 
more negative values of (and therefore faster rotation) 
than the field star sample. Nevertheless, the distribution 
of the sample mean of for the 90% of the SPOCS con- 
trol sample with -5.2 < < 2 is = -0.25 ± 0.19; 
for the host stars of transiting exoplanets an equivalent 
calculation yields = —0.49 ± 0.19. I plot as a func- 
tion of stellar mass for both samples in Figure |4j There 
is a very sharp cutoff in the control sample at » 2 
at all masses, and this sharp cutoff is due to lack of a 
constraint on the inclination of individual stars in the 
control sample. As a result, cr^ im in Equation [2] can be 
large, resulting in small values of 0. 

I find that HAT-P-7, HAT-P-14, HAT-P-16, HD 17156, 
Kepler-5, Kepler-7, TrES-4, WASP-1, WASP-12, and 



WASP-14 have > 2.9, indicating vsmi b s too small to 
be explained by the uncertainties in their masses, radii, 
or ages if each star's spin and its exoplanet's orbit have 
the same inclination to the line of sight. Unless the host 
stars of these transiting exoplanets are slower rotators 
relative to the model prediction than any of the 866 stars 
in the control sample from the SPOCS catalog, the most 
natural explanation for these discrepancies is that the in- 
clination of the spin of these host stars i s is significantly 
different that the inclination of their exoplanet's orbit i p . 
If these anomalously slow projected rotation values are 
the result of spin-orbit misalignment, then the difference 
between v sin i i> s and v sin i sim can be used to compute 
the four possible degenerate values of i s (because sini is 
not one-to-one for 0° < i < 360°). I report these val- 
ues for the potentially spin-orbit misaligned systems in 
Tabled I also identify CoRoT-2, HAT-P-2, HD 189733, 
Kepler-8, and XO-3 as fast rotators. Note that CoRoT-2 
(jAlonso et al.|[20fjl and HD 189733 (pouchy et al.|[200l 
are known to be more active than the Sun, and therefore 
potential fast rotators. HAT-P-2 has been observed to 
have a large radial velocity jitter, which may indicate 
an act ive and possibly quickly rotating star (jBakos et al.l 
[2007a|). 

All ten candidate misaligned systems have host stellar 
mass M* > 1.2 M Q . Again, either these stars have rota- 
tion properties unlike any of the 276 control stars in the 
mass range 1.2 Af Q < M* < 1.5 M Q , or these stars are 
observed pole-on and therefore spin-orbit misaligned. 

2.4. Verification of the Robustness of This Result 

In order to verify the robustness of this result, I have 
repeated my calculations with several different sets of in- 
put parameters. Indeed, there are two potentially uncer- 
tain inputs to my simple empirical model: the power-law 
exponent a in the relation cx t a and the ages of the 
stars in the SPOCS control sample. It is certain that the 
rotation of isolated Sun-like stars slows down with time, 
so a > always. To determine if my result is sensitive to 
a, I have repeated my calculation for a — 1/3, a = 2/3, 
and a — 1. If a = 1/3, then all of the systems identi- 
fied with anomalously slow projected rotation in §2.3 are 
again identified as anomalously slow projected rotators. 
If a = 2/3, eight of the original ten systems are identified 
as anomalously slow projected rotators (HAT-P-7, HAT- 
P-14, Kepler-5, Kepler-7, HD 17156, TrES-4, WASP-12, 
and WASP-14). If a = 1, only three of the original ten 
systems are identified as anomalously slow projected ro- 
tators (HAT-P-7, HAT-P-14, and WASP-14). ' 

I have also redone my calculation using a completely 
uninformative prior on the ages of the stars in SPOCS 
sample. That is, I assumed that the age of each star 
was uniformly distributed between 650 Myr and its main 
sequence lifetime. Again, all ten systems identified in 
§2.3 as anomalously slow projected rotators and therefore 
possibly spin-orbit misaligned are again identified in this 
case. 

Angular momentum exchange between a host star and 
its close-in planetary companion is also unlikely to signif- 
icantly slow the rotation of the host star and produce an 
anomalously slow projected rotation velocity. For that to 
occur, the transfer of angular momentum from host star 
to exoplanet via tidal forces requires the rotation period 
of the star to be shorter than the orbital period of the 
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planet. Indeed, the median orbital period of the transit- 
ing planets orbiting the stars identified as anomalously 
slow projected rotators is P w 3.2 days, while the median 
expected rotation period of their host stars is P* ~ 6.3 
days. 

The presence of slightly evolved stars among the host 
stars of transiting exoplanets will also not produce false 
identifications as anomalously slow projected rotators. 
To see why, note that there are also slightly evolved stars 
in the control sample of SPOCS stars and recall that the 
threshold O required to identify a star as an anomalously 
slow projected rotator was set by the maximum disagree- 
ment between model and observation in the control sam- 
ple. In other words, the presence of slightly evolved stars 
has already been accounted for in setting the threshold. 
I plot in Figure [5] the distribution of the SPOCS sample 
in a theoretical HR diagram along with those transit- 
ing exoplanet hosts with trigonometric parallaxes. I also 
indicate the locations of those stellar hosts that I iden- 
tified as anomalously slow projected rotators and there- 
fore potentially spin-orbit misaligned systems. Figure [5] 
shows that many of the systems I identified as anoma- 
lously slow projected rotators appear to be evolving off 
of the zero-age main sequence (ZAMS). Nevertheless, the 
fact that there are many stars in the control sample of 
SPOCS stars in similar evolutionary states suggests that 
the effects of stellar evolution on the rotation properties 
of evolved stars have already been accounted for in my 
analysis. Moreover, there are several hosts of transit- 
ing exoplanets that have evolved off of the ZAMS that I 
do not identify as anomalously slow projected rotators. 
Collectively, these facts suggest that effects of stellar evo- 
lution cannot explain the rotation properties of the ten 
systems I identified as anomalously slow projected rota- 
tors. 

Additional selection effects in either the SPOCS con- 
trol sample or the transiting exoplanet sample are un- 
likely to affect this result. Indeed, the fact that SPOCS 
stars were preferentially selected to be those stars with 
small enough v sin i to enable high precision Doppler ra- 
dial velocity measurements biases my control sample to 
smaller usini at fixed stellar mass and age. However, 
since the candidate spin-orbit misaligned stars appear to 
have anomalously slow projected rotation, the bias acts 
against the false identification of a star as an anomalously 
slow projected rotator. Though magnitude-limited tran- 
sit surveys are biased toward stars slightly more massive 
than the Sun, there are 276 stars in the control sam- 
ple with > 1.2 M© to compare with. In addition, 
though targeted transit surveys are biased toward less 
active stars, the control sample of SPOCS stars in sim- 
ilarly biased. Transit surveys should also not be biased 
toward evolved stars, as though both transit probability 
and transit duration increase linearly with host stellar ra- 
dius (favoring evolved systems), transit depth decreases 
quadratically with host stellar radius (favoring unevolved 
systems). For that reason, the chance of detecting a tran- 
siting planet is roughly insensitive to host stellar radius. 

2.5. Comparison with Rossiter- McLaughlin 
Measurements 

For the anomalously slow projected rotators identified 
in §2.3, I list in Table [2] the inclination to the line of 
sight of the stellar rotation axis if the anomalous pro- 



jected rotation is due to spin-orbit misalignment along 
the line of sight. I also note published RM measure- 
ments for four of the systems I identify as apparent slow 
projected rotators; HAT-P-7 and WASP- 1 4 has been 
identified as a spin-orbit m isaligned system (|Winn et al.1 
120091* I Johnson et al.ll200l . while HD 17156 and TrES-4 
are cons istent with spin-orbit ali gnment in the plane of 
the sky (|Narita et al.|[2009l l2010h . 

Though HD 17156 is consistent with spin-orbit align- 
ment in the plane of the sky from RM measurements, 
my technique measures misalignment along the line of 
sight. There is other circumstantial evid ence that HP 
17156 is misaligned along the line of sight. Fische r~et al.l 
(2007) measured a usini = 2.6 km s , a stellar radius 
P* = 1.47 Rq and used Ca II H and K to infer a rotation 
period P» = 12.8 days. If the system has i = 90°, then 
a period P* = 12.8 days and the inferred stellar radius 
P* = 1.47 Rq implies vsmi = 5.8 km s _1 (as compared 
to vsmi — 6.9 km s _1 expected based on my simple 
model). Clearly, it's possible that in this case the ac- 
tivity indicator gives a spurious period estimate; still, if 
the discrepancy is due to misalignment then i s = 26° (as 
compared to i s = 22° degrees from my measurement). 

TrES-4 is consistent with alignment in the plane of 
the sky from measurements of the RM effect. However, 
that is a different angle than that measured through my 
analysis. It could be that the geometry of the TrES-4 
system is apparently aligned in the plane of the sky but 
misaligned along the line of sight; alternatively, TrES-4 
could simply be an anomalously slow rotator. 

I find that ten of the 75 systems in Table[T]have anoma- 
lously slow projected rotation and are likely spin-orbit 
misaligned along the line of sight. In contrast, the latest 
Rossiter-McLaughlin measurements reported and sum- 
marized in Triaud et al. (2010, submitted) suggest that 
eight of the 26 systems with RM measurements have sig- 
nificant spin-orbit misalignment in the plane of the sky. 
The higher fraction of systems observed to be misaligned 
by the RM effect is expected, as the RM effect can detect 
much smaller degrees of misalignment than the coarser 
technique described here. Assume for the moment that 
the smallest degree of misalignment this technique can 
identify is \i p — i s \ m in ~ 50° (as suggested by the min- 
imum difference from Table [2]) . Of the 26 systems with 
measurements of the RM effect, six have |A| > 50°. At 
the same time, the slightly higher incidence of misaligned 
systems in the sample of RM measurements might be re- 
lated to the fact that RM measurements are resource 
intensive. That is, systems that seem to have a higher a 
priori probability of spin-orbit misalignment based on the 
current understanding of spin-orbit misalignment (e.g. 
eccentric systems) are more likely to be targeted for RM 
measurement. The fact that this analysis includes all 
known transiting exoplanet systems without the target- 
ing biases inherent in the current sample of RM mea- 
surements may even produce a less biased spin-orbit mis- 
alignment distribution. 

3. DISCUSSION 

I identify ten host stars of transiting exoplanet systems 
in which either the host star is more slowly rotating rel- 
ative to the simple model than any of the stars in the 
control sample of SPOCS stars or in which the star is 
viewed pole-on and the system is spin-orbit misaligned. 
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This degeneracy can be broken for at least Kepler-5 and 
Kepler-7, as high precision Kepler photometry should re- 
veal their rotation periods. If the periods are indeed 
short, then the small observed vsini indicates spin-orbit 
misalignment along the line of sight. Alternatively, if the 
observed rotation periods are long, then Kepler-5 and 
Kepler-7 are systems in which the host star has some- 
how lost angular momentum much more quickly than 
expected by models of Sun-like stellar spin-down. High 
precision photometry from other sources might enable 
similar measurements for the other systems identified as 
anomalously slow projected rotators. 

All ten of the candidate misaligned systems have host 
stars more massive than the Sun. To determine the sig- 
nificance of this observation, I performed a Monte Carlo 
simulation in which I randomly selected with replace- 
ment ten stars from the 75 listed in Table [1] I repeated 
this process 10 6 times. Less than 0.01% of the Monte 
Carlo trials produced a sample in which all ten stars had 
1.2 M© < M* < 1.5 A/©. Therefore, the probability 
that all ten candidate misaligned systems would be iden- 
tified around massive stars by chance is less than one in 
ten thousand. Similarly, the median mass of the planets 
in the candidate misaligned systems is 2 ± 0.6 Jupiter- 
masses, while the median mass of the whole sample of 
75 systems is 1 ± 0.1 Jupiter-masses. Interestingly, there 
may be a hint of the same apparent overrepresentation 
of massive planets in spin-orbit misa ligned systems from 
Rossiter-McLaughlin measurements ([Torres et al.ll2010h . 

Eccentricity also seems to be related to misalign- 
ment probability. Indeed, of the nine transiting systems 
with FGK stellar hosts and eccentricity e > 0.1 (HD 
80606, HD 17156, HAT-P-2, WASP-8, XO-3, HAT-P-11, 
WASP-17, CoRoT-9, and HAT-P-14), there is now evi- 
dence from either measurements of the RM effect or this 
analysis that six of those systems are significantly mis- 
aligned (HD 80606, HD 17156, WASP-8, XO-3, WASP- 
17, and HAT-P-14). The other three systems are unlikely 
to be observed to be spin-orbit misaligned for a number 
of reasons: HAT-P-2 is likely an intrinsically fast-rotator 
and therefore unlikely to be identified as misaligned along 
the line of sight by the method described in this paper 
even if it does have significant spin-orbit misalignment 
along the line of sight, HAT-P-11 may have too small a 
radius for measurements of the RM effect, and CoRoT-9 
has the largest pericenter distance by far of the known 
transiting planets and almost certainly has had a differ- 
ent evolutionary history than the other known transiting 
planets. To determine the significance of the apparent 
correlation between eccentricity and spin-orbit misalign- 
ment, I performed a Monte Carlo simulation to deter- 
mine the expected number of systems with eccentricity 
e > 0.1 from a sample of 16 transiting systems - the to- 
tal number of significantly misaligned systems from RM 
measurements and this analysis - expected if spin-orbit 
misalignment is unrelated to eccentricity. I randomly se- 
lected with replacement the eccentricities of 16 planets 
from the 75 systems listed in Table [TJ and I repeated this 
process 10 6 times. I find that the mean number of sys- 
tems with e > 0.1 in a sample of 16 systems under the 
assumption that eccentricity is unrelated to spin-orbit 
misalignment is 2 ± 2. The observation of six systems 
with e > 0.1 from a random sample of 16 systems occurs 
in less than one in 150 trials, indicating that eccentricity 



and spin-orbit misalignment are related. 

In an attempt to derive some constraint on the under- 
lying spin-orbit misalignment distribution in the 75 sys- 
tems in Table [U I performed two completeness calcula- 
tions for two different i s distributions resulting from two 
different ip distributio ns: (1) corresponding to planet- 
planet scattering fe.g. iChatteriee et al.l 12008). a case in 
which the inclination of the host star is completely in- 
dependent of the alignment of the orbit of its planet 
i s ~ arccos(l — U), where U ~ Unif (0, 1); and (2) 
corresponding to Kozai cy cles and tidal friction (e.g. 
iFabrvckv fc Tremaind 1200^ 1 . a case in which the distri- 
bution of i s is derived from the ijj dis t ributi on as given 
in Figure 10 of IFabrvckv fc Tremaind" (|2007f ) and the A 
distribution as given in the RM compilation in Table 5 of 
Triaud et al. (2010, submitted). Assuming a threshold 
detectable misalignment that corresponds to the smallest 
degree of misalignment identified (\i p — i s \ > 50°) in 75 
systems, I find that the expected number of detections 
in case (1) is 17 ± 4 and in case (2) it is 27 ± 4. Even 
though each individual system in this analysis is securely 
identified as an anomalously slow projected rotator and 
therefore possibly spin-orbit misaligned, both complete- 
ness calculations are strongly dependent on the smallest 
degree of misalignment this technique can identify. That 
quantity is uncertain, and as such more data is neces- 
sary before any conclusive statements about distribution 
of spin-orbit alignments is made. The full sample of Ke- 
pler transiting planet detections and host star properties 
available at the end of its four year mission might resolve 
this issue. 

Collectively, this analysis indicates that predictions 
of spin-orbit misalignment for a population of close-in 
planets entirely produced by planet-planet scattering or 
Kozai cycles with tidal friction overpredict the number 
of misaligned systems. As a result, there seems to be 
two populations of close-in planets: a population that 
is spin-orbit aligned and a population that is apparently 
spin-orbit misaligned. The processes that lead to mis- 
alignment or to the survival of misaligned systems seem 
to operates most efficiently in systems with massive host 
stars and planets. Indeed, the striking appearance of can- 
didate misaligned systems at M* > 1.2 M in Figured] 
indicates that the processes that lead to spin-orbit mis- 
alignment and survival are threshold processes, similar to 
the rapid increase i n exoplanet incidence with host stel- 
lar m etallicity (e.g. Santos et al. 2004; Fische r fc Valentil 
2005). In this case, the transition to frequent appar- 
ent spin-orbit misalignment occurs at the same stellar 
mass at which Sun-like stars with near solar metallic- 
ity develop radiative envelopes. This dramatic change 
in stellar structure will likely strongly effect the poorly- 
understood tidal processes that play a major role in the 
formation, evolution, and long-term survival of spin-orbit 
misaligned systems. 

4. CONCLUSION 

There are ten transiting exoplanet systems from a 
larger sample of 75 systems in which the projected rota- 
tion velocity of the host star assuming spin-orbit align- 
ment is apparently smaller relative to a simple empirical 
model of stellar spin-down than any of 866 stars in a 
representative control sample of stars in the same mass 
range from the Spectroscopic Properties of Cool Stars 
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(SPOCS) catalog: HAT-P-7, HAT-P-14, HAT-P-16, HD 
17156, Kepler-5, Kepler-7, TrES-4, WASP-1, WASP-12, 
and WASP-14. Unless the host stars of these transit- 
ing exoplanets are slower rotators relative to the model 
prediction than any of the 866 stars in the representa- 
tive control sample, the anomalously low v sin i values 
are best explained by a pole-on view of the star and 
therefore spin-orbit misalignment. I only identify signif- 
icant possible spin-orbit misalignment in systems with 
host stars more massive than the Sun, as in all ten sys- 
tems the host stars have stellar masses M* in the range 
1.2 Mq < M* < 1.5 Mq. This observation is significant, 
as there is less than a one in ten thousand chance of its oc- 
currence by chance. There is evidence that systems with 
massive and eccentric planets are also those most likely 
to be significantly spin-orbit misaligned. Given current 
measurements of the Rossiter-McLaughlin effect and fre- 
quency of misalignment along the line of sight determined 
in this analysis, models of close-in planet formation in- 
voking planet-planet scattering or Kozai cycles and tidal 
friction alone overpredict the number of misaligned sys- 
tems. As a result, there seems to be two populations 
of close-in planets: those that are aligned and those that 
are apparently misaligned. The process that leads to mis- 
alignment or survival of misaligned systems likely oper- 
ates more efficiently in systems with both massive stars 
and planets. Interestingly, apparently misaligned sys- 
tems appear to occur preferentially in systems in which 
the envelope of the host star is radiative, and this tran- 
sition in stellar structure is likely related through tidal 
processes to the formation, evolution, and long-term sur- 
vival of spin-orbit misaligned systems. 

I thank Greg Laughlin and Connie Rockosi for useful 
comments and conversation. This research has made use 
of NASA's Astrophysics Data System Bibliographic Ser- 
vices. This material is based upon work supported un- 
der a National Science Foundation Graduate Research 
Fellowship. 
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Figure 1. The geometry of transiting exoplanet systems following the conventions described in Fabrycky & Winn (2009), with the z-axis 
oriented toward the observer. The exoplanet orbit is defined by the vector n p and is inclined at an angle i v to the observer's line of sight. 
The host stellar spin is defined by the vector n s and is inclined at an angle i s to the observer's line of sight. Transit measurements can 
determi ne i p , and measurements of the Rossiter-McLaugh lin effect can reveal A, the angle in the xy-plane between the projection of n p 
and n s (Rossiter 1924; McLaughlin 1924; Ohta ct al. 2005). The angle i s is difficult to measure for individual stars, but as I show in this 
paper, it may be possible to identify systems in which i s and i v are misaligned by a large degree. 
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Figure 2. Stellar rotation period of Sun-like stars as a function of mass and ag e. Though a population of Sun-like stars is formed with a 
wide range of rotation periods P» ( Attridgc & Hcrbst 1992; Choi & Hcrbst 1996), by the time that population is the age of the Hyades and 
Praesepe period is only a function of mass (e.g. Tassoul 2000; Irwin & Bouvier 200!^), Main sequence stars lose angular momentum an d 
evolve from from this initial condition, likely through a mag netized stellar wind (e.g. [Weber & Davis 1967; Mcstcl 1968; Kawalerl 119881) . 
As a result, stellar periods increase on the main sequence as they age (e.g. Kraft 1967) as P„ <x t 1 / 2 , in accordance with the empirical 
relation vsvai oc t -1 / 2 note d by Skumanicb| Q1972I). I binned the avail able Hyades and Praesepe data ijRadick et al. 1987; Pross er" et al.l 
11991 ISchoIz & Ei sloffcl 2007]) as summarized in Irwin & Bouvier (2009) (black points) in bins of 0.1 Mq (blue squares), supplemented it 
with average rotation values for higher mass field stars ((McNally 1965) (orange squares), and used natural cubic spline interpolation to 
compute a relationship between mass and period for stars between 0.2 Mq < M„ < 1.75 Mq (solid gray curve). From this initial condition, 
I evolve the rotation period as P, oc t 1 / 2 (dotted gray vertical lines). The length of the vertical lines correspond to the main sequence 
lifetime of a star of the given mass, and I denote Gyr intervals (from 1 Gyr to 13 Gyr) with gray triangles along each vertical line. 
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Figure 3. The distribution of the rotation statistic for stars that host transiting exoplanets (solid blue line) compared to a control 

sample from the SPOCS catalog of I Valenti &; Fischerl 120051) (solid gray line). I define = («sini s i m — v sin i i, 3 ) / y + a \i m such 

that stars that have a vsini smaller (larger) than expected given their mass, age, radius, and the uncertainties in those parameters have 
positive (negative) values of 0. By comparison with the control sample, transiting exoplanet host stars that have > 2.9 have vsini 
so small that they cannot be explained by evolving the Hyades and Praesepe data according to the P* oc t 1 / 2 relation given the system 
parameters and uncertainties and assuming alignment between the exoplanet orbit and the spin of its host. I argue that unless these stellar 
hosts of transiting exoplanets are slower rotators relative to the model than any of the 866 stars in the SPOCS control sample, the most 
natural explanation for this result is that i s and i p are misaligned. 
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Figure 4. The rotation statistic © as a function of stellar mass for both stars that host transiting exoplanets (blue points) and control 
stars from the SPOCS catalog (gray points). Note that stars with © > 2.9 - anomalously slow projected rotation assuming spin-orbit 
alignment - are preferentially massive. In the sample of 75 systems in Table [T] the probability that all ten misaligned systems have host 
stellar mass in the range 1.2 Mq < M* < 1.5 Mq by chance is less than one in ten thousand. In addition, the exoplanet host stars with 
> 2.9 host planets with a range of semimajor axes, not just those at small semimajor axes where tidal forces are most important. In any 
case, the rotation period of each exoplanet host star with > 2.9 is expected to be longer than the orbital period of its planet, so those 
stars are unlikely to transfer angular momentum to their planets. 
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Figure 5. A theoretical HR diagram of the control sample of SPOCS stars (gray points) overplotted with the anomalously slow projected 
rotators identified in this analysis. I also plot other hosts of transiting cxopl anets with pr ecise trigonometric parallaxes (black points). 
The dashed lines are Padova isochrones for ages 0.1, 3, 5, and 10 Gyr {Marigo ct al. 200^). There are both young and evolved stars in 
the control sample, so the SPOCS sample is a fair control sample in regards to evolutionary state. In addition, the fact that there are 
evolved hosts of transiting exoplanets that are not identified as anomalously slow projected rotators suggests that the rotational evolution 
of evolved main sequence stars alone cannot explain the anomalously slow projected rotation identified in some systems. 
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Table 1 

Transiting Exoplanct System Properties 
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0.50 
1.00 
1.00 
1.00 
1.00 



0.58 
0.30 
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0.50 
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1.00 
0.50 
0.50 
1.00 
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0.40 
1.00 
0.50 
0.50 
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0.50 
0.05 
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0.70 
1.00 
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1.00 
1.00 
0.50 



1.50 
1.00 
1.00 
1.00 
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1.50 

0.57 
2.00 
1.00 
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0.60 
0.60 
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1.60 
0.60 
10.00 
0.67 
0.30 
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Table 1 — Continued 



System 


M, 






TV 


vsini oba a oba 




V Sin l s i m °~sim 





Reference 




[Me] 


[Re] 


[Gyr] 


[Gyr] 


[km s _1 ] [km s — 1 ] 


[deg] 


[km s _1 ] [km s — 1 ] 






XO-5 


1.00 


l.n 


7.70 


9.30 


0.70 0.50 


86.80 


1.71 0.05 


2.01 


87,88 



Note. — (DIBarge et alj POOSTl; (2)IAlonso et alj POOST) : f3)[BouchYeF al. (2008); (4) Dclcuil ct al. (2008); (5) Aigrain ct al. (200 jg); (6) 
IMoutou etldJ1200Sll; ( 7) IRauer et alj 12009 ) : (8) Fridlund et al.l (120101); (9)ILeger etaLl (12009 1); (l O)IDeeg et al J 12010); (ll) ISantos et al. l 
(2004); (12) Torres ct al. (2008); (13) Bean et al. (2008); (14) Charbonncau ct al. (2009); (15) Johnson ct al. (2008); (16) Pal ct al. 
(2010); (17) Torres ct al. (2007); (18) Kovacs ct al. (2007); (19) Bakos ct al. (2007b); (20) Novcs ct al. (2008); (21) Pal ct al. (2008); 
(22)ILatham et a l. (2009); (23) Shpo rer et alj (120091) ; (24)IBakos et alj d20l6Tj ~ (25) Hartman ct al. (2009); (26) Bakos ct al. (2009% (27) 
[ Torres et all 12010) ; (28) Buch have"et alj ( 120101 ) ; (29) [Fischer et alj (1 200711; (30) I Winn et alj 1 2009a); (31 )IFischer fc Val cnti (2005 ); (32) 
Pont ct al. (2009); (33) Fosscy ct al. (2009); (34) Sato ct al. (2005); (35) Nutzman ct al. (2009); (36) Winn ct al. (2007); (37) Triaud ct al. 
( 12009); (38)IWinn et al. (2005); (39) Borucki ct al. (2010); (40) Koch ct al. (2010); (41) Dunham ct al. (2010); (42) Latham et alj d20ld) T 
(43) IJenkins et alj 112010) ; (44) IMelo et alj 12009) ; (45) ISouthworthl d2008 1)~(4 6) [Hillon et alj 12000) ; (46) ISnellen et alj 12009) ; (48) 
Laughlin ct al. (2005); (49) Narita ct al. (2007); (47) Sozzctti ct al. (2007); (51) Daemon ct al. (2009); (52) Sozzctti ct al. (2009); (53) 
Cameron ct al. (2007); (54) Pollacco ct al. (2008); (55) Gibson ct al. (2008); (56) Gillon ct al. (2009a); (57) Southworth ct al. (2009a); 
(58)IGilIon et alj 12009bl) ; (59) Hcllicr ct al. (2009a); (60) Quelo z et al. (201 , subm itted ); (61) IChristian et al.l 12009); (62) Bako s et all 
(2009a); (63) Hcbb ct al. (2009); (64) Skillcn ct al. (2009); (65) Joshi ct al. (2009); (66) Johnson ct al. (2009); (67) West ct al. (2009); 
(68) Lister ct al. (2009); (69) Anderson ct al. (2010);' (70) Hcllicr ct al. (2009b); (71) Southworth ct al. (2009b); (72) Hcbb ct al. (2010); 
(73) Bouchy et al. (2010, submitted); (74) Maxted et al. (2010, submitted); (75) Street et al. (2010, submitted); (76) Enoch et al. (2010, 
submitted); (77) Smalley et al. (2010, submitted); (78) W est et al. (2010, submitte d); (7 9) Hellier et al. (2010, s ubmitt ed); (80) Cameron 
et al (2010 , submitted); (81) [McCullough ct al. (2006); (82) Bur ke et aLl 120071): C83 ) [Fernandez ct al. (2003); (84) [Jolms-Krul l et al] 
120081) ; (85) IWinn eTaU 120081) ; (86) IMcCullough et al.l pOOSl) ; (87) IBurke et alj 120081) ; (88) IPal et al.l 12009b!) 
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Table 2 

Misalignment Properties 



oystem 


%p 
[degj 


is,i 
[degj 


is,2 
[degj 


*s,3 
[degj 


*s,4 
[deg] 


A 

[deg] 


Reference 


HAT-P-7 b 


85.70 


6.7 


170 


190 


350 


-177.5 


Winn et al. (2009b) 


riAl-F-14 D 


83.50 


21 


160 


200 


340 




HAT-P-16 b 


86.60 


34 


150 


210 


330 






HD 17156 b 


86.20 


22 


160 


200 


340 


10 


Narita et al. f2009) 


Kepler-5 b 


86.30 


19 


160 


200 


340 




Kepler-7 b 


86.50 


20 


160 


200 


340 






TrES-4 


82.86 


31 


150 


210 


330 


-6.3 


Narita et al. (2010) 


WASP-1 b 


88.25 


36 


140 


220 


320 




WASP-12 b 


83.10 


8.8 


170 


190 


350 






WASP-14 b 


84.32 


15 


160 


200 


3 10 


-33.1 


Johnson et al. ("20091 



Note. — Column 2 is the inclination of the transiting planet's orbit to the line of sight i p . Columns 3 through 6 are the four possible 
values of i s given the value of sin i inferred from the difference in tisini {, s and vs'mi s i m . Column 7 is the measured value of A from the 
Rossiter-McLaughlin measurement taken from the indicated reference. 



